Translations and Rotations are correlated in Granular Gases 
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In a granular gas of rough particles the axis of rotation is shown to be correlated with the 
translational velocity of the particles. The average relative orientation of angular and linear velocities 
depends on the parameters which characterise the dissipative nature of the collision. We derive a 
simple theory for these correlations and validate it with numerical simulations for a wide range 
of coefficients of normal and tangential restitution. The limit of smooth spheres is shown to be 
singular: even an arbitrarily small roughness of the particles gives rise to orientational correlations. 



PACS numbers: 45.70.-n,45.70.Qj,47.20.-k 

Dilute systems of macroscopic particles, called Gran- 
ular Gases, show many novel and surprising phenom- 
ena when compared to molecular gases. The particles 
of Granular Gases are macroscopic bodies which in gen- 
eral dissipate energy upon collision. As a consequence, 
such gases demonstrate features, which drastically dif- 
fer from molecular gases: The velocities are not dis- 
tributed according to a Maxwell-Boltzmann distribution 
[1 II H II II II Ql; equipartition does not ho ld 11 and a 
homogeneous state is in general unstable 0, lid lllj . In 
this Letter we present another unexpected result: 

The angular and linear velocities of rough particles 
are correlated in direction. In dependence on the coeffi- 
cients of restitution, characterizing the dissipative parti- 
cle properties, the rotational motion may be preferably 
perpendicular to the direction of linear motion, similar 
to a sliced (spinning) tennis ball, or in parallel to it, like 
a rifled bullet. Fig. ^ Surprisingly, the limit of vanishing 
dissipation of the rotational motion does not exist, that 
is, any arbitrarily small roughness leads to a macroscopic 
correlation between spin and velocity. We present a ki- 
netic theory that quantifies this new effect and find good 
agreement with large scale numerical simulations. It is 
expected that the reported correlation between spin and 
linear velocity may have important consequences in un- 
derstanding natural Granular Gases, such as dust clouds 
or planetary rings. 

Model- We consider a monodisperse granular gas con- 
sisting of N hard spheres of radius a, mass m and moment 
of inertia / = qma^ {q — 2/5 for homogeneous spheres 
addressed here) . The degrees of freedom are the particles' 
position vectors {vi}, translational velocities {v.^} and ro- 
tational velocities {w;} for i = 1,2,...N. The dynamic 
evolution of the system is governed by instantaneous two 
particle collisions, such that in general both translational 
and rotational energy is dissipated. Introducing the rel- 
ative velocity at the point of contact 

g = vi -■V2 + an X {lji + LJi) (1) 



the collision rules specify the change of g in the direction 
of n = (ri — r2)/|ri — and perpendicular to n: 

(g • n)' = -e„(g • n) (g x n)' = +et(g x n), (2) 

where the primed values refer to the post-collision quan- 
tities. The coefficients of restitution in normal and tan- 
gential direction, < e„ < 1 and — 1 < e* < 1, charac- 
terize the loss of energy and, thus, describe the slowing 
down of the linear and rotational motion of the parti- 
cles. These coefficients are the central quantities in the 
Kinetic Theory of granular gases 12] . For e„ = 1 (elastic 
spheres) and et = ±1 (perfectly smooth/rough spheres) 
the energy is conserved, while for e„ — > and — > dis- 
sipation is maximal. Combining the collision rules (Eq. 
(P|)) with conservation of angular and linear momentum, 
allows one to express the post-collision velocities in terms 




FIG. 1: In granular gases, depending on the coefficients of 
restitution, spin (red arrow) and linear velocity (grey arrow) 
are oriented either preferably in parallel (cannon ball) or per- 
pendicular (sliced tennis ball). This is a fundamental differ- 
ence to molecular gases, where spin and velocity are uncorre- 
lated in orientation. 
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of the precoUisional ones 

v| — vi — d , Lj{ — LJi + ( —] n X d, 

\aqj 

■Vn = V2+S, o)^ = a;2 + 1 — I n X 5, (3) 

\aqj 

with 6 = r]tg + (r/„ - r7f)(n ■ g)n and ?]„ = (1 + e„)/2, 
rjt^qil-et)/2{l + q). 

When particles coUide according to the above laws, two 
processes take place: (i) dissipation of energy and (ii) ex- 
change of energy between the rotational and translational 
degrees of freedom. The first process is described by two 
time-dependent granular temperatures 

N I ^ 

^- = ^E^' T„, = — (4) 

i—l i—1 

The second process drives the system to a quasi steady- 
state, which is characterized by a constant ratio of the 
two temperatures, r — Trot/Ttr — const. In this state 
both temperatures continue to decay with their rates tied 
together by Trot/Ttr = r. The ratio of temperatures, r, 
depends on the coefficients of restitution e„ and et and 
can take values smaller or larger than one. If r < 1 
collisions damp the translational motion more efficiently 
than the rotational one, \Ttr\ = iTrotl/f > \Trot\, whereas 
for r > 1 the rotations are damped more efficiently. 

At first glance there is no reason to expect that linear 
and angular velocities v and u) are correlated. Neverthe- 
less, the exchange of energy between the rotational and 
translational degrees of freedom may build up such cor- 
relations. We quantify them by the mean square cosine 
of the angle between v and a;; 

If there are no correlations between the rotational and 
translational motion, as in molecular gases, one obtains 
K = 1/3; in granular gases we find that K in general 
deviates significantly from 1/3. 

The full quantitative understandig of these correlations 
requires detailed mathematical analysis of the collision 
rules Eq. © (see below). Nevertheless simple physi- 
cal arguments are helpful for a qualitative understand- 
ing of these correlations. The transfer of energy from 
the translational to the rotational degrees of freedom 
(and vice versa) depends sensitively on the relative ori- 
entation of Vi — Vj and uJi + ojj. To sustain a quasi 
steady state with a fixed ratio, r, of Ttr and Tiot re- 
quires that ffuctuations away from a given r are effec- 
tively suppressed. In the limit of nearly smooth par- 
ticles, £t ^ 1, the argument is particularly simple: In 
this case, we have r — Trot /Ttr ^ 1 0|, that is, for 
most of the collisions the relative velocity of the parti- 
cles at the point of contact is dominated by rotations 



g = Vi — Vj -I- an X {ijJi + ~ an X (w^ + Ljj). For 
nearly smooth particles the rotational velocities change 
only slightly upon collision, u>l w u)i and « a;^, so 
that the collision rule Eq. © simplifies to 

v/ - Vi « -?7tg « --qta (n x a;^ + n x u)'^) , (6) 

where 77^ ~ (1— e^) ^ 1 . The first term on the right-hand 
side gives a contribution to v/ that is always perpendic- 
ular to the angular velocity us'^. The second term (of the 
same order of magnitude), ?ytan x uj'p has no preferred 
orientation with respect to a;^, - the two vectors u;^ and 
(jj'j being uncorrelated in a dilute gas. henec, the sum of 
all contributions leads to being preferably perpendic- 
ular to 

To describe these phenomena quantitatively beyond 
the limit of nearly smooth spheres, we develop an an- 
alytical theory which is based on an ansatz for the iV- 
particle distribution function. We assume homogeneity, 
except for the excluded volume interaction, and molecu- 
lar chaos, implying that the A^-particle velocity distribu- 
tion factorizes into a product of one particle distributions 

pi(v,a;,t) ^ e"5ivw e^^^-tO [l + 6(t)v^w^F9(cos 6)] . 

(7) 

This ansatz takes into account for the first time orienta- 
tional correlations between linear and angular velocities. 
To lowest order these correlations can be characterized 
by the second Legendre polymonial ^2(2;) = fa;^ — \- 
In general there will be higher-order terms in (cosO), as 
well as non-Gaussian corrections for the velocity and an- 
gular velocity distribution. Here we ignore these terms 
since the simplest ansatz already captures, even quanti- 
tatively, the correlations of interest. The ansatz contains 
three functions, Ttr{t), Trot(i) and &(i), which have been 
calculated self-consistently using the Pseudo Liouville op- 
erator. These functions determine the correlation of Eq. 
(0) according to 

Results- We find that very strong dynamic correlations 
may develop, depending on the initial conditions. Af- 
ter a transient period, the correlation factor K reaches 
a steady-state value which is independent of the initial 
conditions. In Fig. |2] we show this steady-state value 
as a function of the coefficient of tangential restitution 
£t for several values of the normal coefficient £„. We 
compare results from the analytical theory with data 
obtained with Direct Simulation Monte Carlo (DSMC) 
[1^ of = 2 X 10^ particles. The analytical results are 
in good agreement with simulations, especially for small 
dissipation. As can be seen in the figure, the most pro- 
nounced correlations are observed for Et — 0, that is, for 
maximal damping of the tangential motion. The correla- 
tions weaken when the damping decreases with increasing 
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FIG. 2; Steady-state value of the correlation factor K — 
(cos^ 0i) as a function of the coefficient of tangential restitu- 
tion, Et, for different values of normal restitution, e„. Lines 
- analytical theory, points - DSMC simulations. The isolated 
point at the right border where et = 1 indicates vanishing 
correlations for ideally smooth particles. 



et| and are small when the dissipation of the tangential 
and normal motion is comparable, i.e. e„ ~ et- With 
still increasing magnitude of et, the correlations increase 
again and K deviates significantly from the uncorrelated 
value 1/3. 

The complete dependence of the correlation factor K 
on both coefficients of restitution is shown in Fig. El 
along with the contour plots of the temperature ratio 
r. In agreement with the qualitative argument discussed 
above, we see that translational and rotational veloci- 
ties are preferentially prependicular {K < 1/3) in those 
regions where r strongly deviates from 1. 

The right-hand side of Fig. [3 reveals that the limit 
of smooth spheres, is not continuous - contary to naive 
expectations. For perfectly smooth particles, et = 1, 
we obtain K = 1/3, as expected for a molecular gas. 
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FIG. 3: Analytical predictions for the steady-state value of K 
as a function of e„ and et- The black lines show the manifold 
K = 1/3 (no correlations), the white lines are the contour 
plots for r = Trot/Ttr. 



This happens, because translational and rotational mo- 
tion decouple, so that all spins simply persist in time. 
Hence, trivially, angular and linear velocities are not cor- 
related for £t = 1. Nevertheless, even for vanishingly 
small roughness et ^ 1 the correlation factor K notice- 
ably deviates from 1/3. Our analytic theory predicts 



lim (cos^ 8) 



3 (1 - e„) 
8 (7 - e„) ' 



(9) 



in good agreement with DSMC. In the limit et ^ 1 
the relaxation time to the quasi steady state diverges, 
because the exchange of energy between rotational and 
translational degrees of freedom becomes more and more 
inefficient as their mutual coupling vanishes. 

Fig. ^ gives an illustrative example, how correlations 
develop in time. We have chosen parameter values and 
initial conditions such that both regimes - dominance 
of translational motion, r ^ 1 and dominance of rota- 
tional motion, r 1 - are visible. Initially K — l/Z 
and r — > 0, which means that correlations are lacking 
and the particles' spins are vanishingly small. The evo- 
lution proceeds in three stages. In the initial stage rota- 
tional motion is generated mainly in grazing collisions so 
that the particles rotate around an axis perpendicular to 
their linear velocity, like a spinning tennis ball, implying 
(cos^(0)) < 1/3. Once the rotational motion has become 
comparable in magnitude to the translational motion, the 
correlations are small. In this intermediate time regime 
the system is still far from the quasi stationary state 
which, for our choice of e„ and et, is characterized by 
a stationary value of r ^ 1 . It is plausible that correla- 
tions are small in this intermediate regime, because both 
translational and rotational velocities contribute about 
equally to the momentum transfer in collision. Finally in 
the asymptotic state the rotational motion is consider- 
ably more intense than the translational motion, so that 
the quasi stationary state with r ^ 1 is characterized by 
significant correlations K < 1/3. 

Methods- Two different numerical methods were used: 
Direct simulation Monte Carlo (DSMC) Q and Event- 
driven Molecular Dynamics (MD) 15] . In contrast to the 
latter method, where the trajectory of each particle of 
the ensemble is directly computed according to the basic 
kinematics and collision rules, the DSMC is based on the 
solution of the kinetic Boltzmann equation. Because the 
gas is explicitly treated as uniform in this method, and 
spatial correlations between grains are ignored, it allows 
to handle up to ~ 10^ particles, which is necessary for 
the precise measurement of the effect of interest. While 
the MD method is more appealing from a physical point 
of view, the DSMC is much more powerful for a homo- 
geneous granular gas; in the limit of low density the two 
methods are, in principle, identical [lflj |. 

In the analytic approach one can use either the Pseudo 
Liouville operator technique (for details see e.g. 0|) 
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FIG. 4: Relaxation of the temperature ratio r and correlation 
factor K to their steady-state values for £„ = 0.8 and £t = 0.8. 
The arrows indicate the relevant axes. Full lines: analytical 
theory, dots: simulations of = 8000 particles. In the initial 
state Trot = 0, Ttr = 1, K = 1/3. Note the three different time 
regimes: First, from r = to r = 1, very strong correlations 
develop; second, for r ~ 1 correlations are small; finally in 
the quasi steady state with r ^ 1, correlations acquire a 
stationary value K < 1/3. The dashed line shows K — 1/3 
as expected for a molecular gas. 



Details of the calculation will be published elsewhere. 

Conclusions- In conclusion, we reveal a novel phe- 
nomenon, which is unique for granular gases and sur- 
prising in two respects: a) Except for very special values 
of the coefRcients of restitution, St), the linear and 
angular velocities are noticeably correlated. For most of 
the parameter values v and u; are preferably perpendic- 
ular and K < 1/3 like for a sliced tennis ball. In a small 
region of low dissipation they are preferably parallel, so 
that K > 1/3 like for a rifled bullet, b) The limit of van- 
ishing tangential dissipation, — > ±1 is not continuous. 

These results have important consequences for the hy- 
drodynamic theory of dilute granular flows. It was re- 
cently shown that the angular velocity needs to be 
included in the set of hydrodynamic fields. In view of the 
singular nature of the limit of vanishing roughness, per- 
turbation expansions around the smooth limit are ques- 
tionable. Orientational correlation between spin and lin- 
ear velocity presumably also affect the stability of the 
sytem towards shear fluctuations which constitute the 
dominant instabilty of granular flows of smooth particles. 



or equivalently the Boltzmann equation to derive self- 
consistent equations for Ttr{t),Trot{t) and h{t). In the 
quasi stationary state the relative orientation is given by: 
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(cos^ e)c 



~ 5 - 40C +{A- C)^ + (40B + BW)r*-^ ' 

Here r* is the stationary ratio of temperatures which is 
given in(0) together with the coefRcients A,B and C. 
The remaining coefficients are explicitly given by: 
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